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Abstract. The FASTSUM collaboration has been carrying out simulations of A/ = 2 + 1 QCD at nonzero temperature in 
the fixed-scale approach using anisotropic lattices. Here we present the status of these studies, including recent results for 
electrical conductivity and charge diffusion, and heavy quarkonium (charm and beauty) physics. 
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INTRODUCTION 

The heavy-ion collision experiments at RHIC and LHC have given a wealth of information about QCD matter at high 
temperatures, and have provided strong evidence for the existence of a deconfined state of matter, the quark-gluon 
plasma (QGP). Furthermore, it has been shown that at the energy densities reached in these experiments, this state of 
matter can be described well by nearly ideal hydrodynamics, and must hence be strongly coupled. A large amount of 
information has also been obtained about various probes of the QGP, including electromagnetic probes, high-energy 
jets and heavy quarks. 

One of the outstanding challenges is to obtain a clear theoretical understanding of all the phenomena observed 
in the heavy-ion collisions. In particular, first-principles calculations of transport coefficients such as viscosity and 
conductivity would place the hydrodynamical description on a solid footing. Heavy quarkonia have been proposed 
as “thermometers” of the QGP, but this requires a precise knowledge of their dissociation rates as a function of 
temperature. 

Lattice QCD is the method of choice for obtaining quantitative theoretical predictions for QCD matter at high 
temperature and zero or small net baryon density. However, to obtain information about real-time quantities such 
as transport coefficients and dissociation rates, analytical continuation of the imaginary-time information obtained 
from lattice simulation is necessary. This is in principle an ill-posed problem, but it can be addressed using Bayesian 
methods such as the maximum entropy method (MEM) [1] and the new bayesian reconstruction method of [2], or with 
alternative, model independent methods [3, 4]. 

All these methods require a fine resolution in the temporal direction to obtain reliable results. This may be done by 
using anisotropic lattices with a smaller spacing in time than in space. The FASTSUM collaboration has for a number 
of years been carrying out simulations on anisotropic lattices with dynamical fermions with this in mind. Here we 
present some of the most recent results from these simulations. 


TABLE 1. Parameters for the 1st and 2nd generation ensembles. ^ = 
is the anisotropy; Ns is the number of spatial sites and Lj is the 
extent of the lattice in the spatial directions. 
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2.94 







32 

3.94 


TABLE 2. Temporal lattice extents Nt and tempera¬ 
tures T used in our simulations. 


Ar 

Gen 1 

T (MeV) 

T/Tc 

Nr 

Gen 2 

T (MeV) 

T/Tc 

80 

92 
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0.84 
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SIMULATION AND ANALYSIS DETAILS 


We have generated two sets of anisotropic lattice ensembles using improved Wilson-type fermions, see Table 1 for 
details. The 1st generation ensemble [5, 6] had Nf = 2 active flavours with nij^ « 500MeV, while the 2nd generation 
ensemble has Nf = 2 + 1 active flavours with a physical strange quark and « 400MeV. The 2nd generation 
parameters are the same as those used by the Hadron Spectrum Collaboration [7, 8], who have kindly allowed us to 
use their zero-temperature configurations for the purpose of these studies. The 2nd generation configurations as well 
as light and charm quark correlators were produced using the Chroma software system [9] with BAGEL optimisation 
[ 10 ]. 

For each parameter set we have generated gauge configurations at a range of temperatures, see table 2 for details. 
Note that the second generation includes temperatures both below and above the deconfinement transition. The 
pseudocritical temperature Tc was determined from the inflection point of the renormalised Polyakov loop, see [11, 12] 
for details. 

Transport properties and spectral information are encoded in spectral functions p((0,p), which are related to the 
euclidean correlators Ge{t,p) that can be computed on the lattice through the integral relation 
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for a system in thermal equilibrium. We have used the maximum entropy method with Bryan’s algorithm and the 
modified kernel proposed in [13] to determine the most likely spectral function given our correlators. We have also 
employed the novel Bayesian method developed in [2], and will present results using this method in the beauty sector. 


CONDUCTIVITY, CHARGE SUSCEPTIBILITY AND CHARGE DIFFUSION 


The electrical conductivity is related to the low-frequency limit of the current-current spectral function, 
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FIGURE 1. Left: light quark vector current spectral functions for three different temperatures below, near and above the 
pseudocritical temperature. The bands represent the statistical uncertainty. Right: the electrical conductivity (including contributions 
from u, d and s quarks) as a function of temperature. The boxes represent the uncertainties due to variations in the default model, 
while the error bars represent the total (systematic and statistical) uncertainty. 


where = L/ /) jfi is the electromagnetic current operator and qf are the quark charges in units of the elementary 
charge e. We have used the exactly conserved vector current on the lattice to compute Gf“(T), see [14, 12] for details. 
In order to compare results obtained from simulations with different flavour content, in particular results from 2 (u, d) 
and 3 (u, d, s) flavours, the spectral functions and conductivity have been divided by the factor Cem = i-®- 

the sum of the squares of the quark charges. 

We have also computed the susceptibilities of baryon number, isospin and electric charge as a function of tempera¬ 
ture [15, 12]. These are of substantial phenomenological interest as they are related to the event-by-event fluctuations 
of these quantities in heavy-ion collisions, and may be used to locate the transition to the QGP and the critical endpoint 
in the temperature - chemical potential plane, if it exists. 

Furthermore, the conductivity is directly related to the charge diffusion coefficient Dq through the relation Dq = 
<yIX q, where Xq is the electric charge susceptibility. While the analogous diffusion coefficient for heavy quarks has 
been widely studied, much less is known about diffusion of electric charge carried by light quarks. Here we present 
the first lattice calculation of this quantity. All the results in this section have been obtained with our 2nd generation 
ensemble. 

In the left hand panel of figure 1 we show the light (u, d) quark contribution to Pem{(o) as a function of co for three 
different temperatures: below, near and above the pseudocritical temperature Tc. Below Tc we can clearly identify the p 
meson peak in the spectral function, while this peak disappears at higher temperatures, in accordance with the p meson 
no longer being bound. Above Tc we find that p{(o)/(0 has a nonzero intercept with the 0 = 0 axis, signalling a nonzero 
value for the conductivity. The MEM analysis has been carried out using a default model m{co) = moco{b + Co), which 
allows for a finite, nonzero transport contribution while reproducing the continuum free spectral function at large co. 

The resulting values for the conductivity are shown in the right hand plot of figure 1. We have studied the stability 
of our results with respect to variations in the default model parameter the resulting systematic uncertainties are 
represented by the filled boxes in the figure, while the error bars represent the total (systematic and statistical) 
uncertainties. We have also investigated the stability of our results with respect to other systematics including the 
range and number of time slices included, and found that they are stable in all cases. 

In the left panel of figure 2 we show our results for the isospin, electric charge and baryon susceptibilities as 
functions of temperature. They all show a similar behaviour, increasing rapidly near the crossover temperature Tc and 
approaching the value for a gas of free quarks and gluons at high temperature, as expected. 

The right panel of figure 2 shows the charge diffusion coefficient Dq as a function of temperature. Our results 
suggest that it has a minimum near Tc, with a value close to that predicted in AdS/CFT models, Dq = 1 /{2tzT). 













T [MeV] T [MeV] 


FIGURE 2. Isospin, electric charge and baryon susceptibility XuXq^Xb, normalised by their values x^^ in the noninteracting 
case (left) and charge diffusion coefficient Dq (right) as function of temperature. These results include the contributions from u, d 
and s quarks. The meaning of the filled boxes and error bars are the same as in figure 1. 




FIGURE 3. Charmonium spectral functions from the second generation ensemble. Pseudoscalar channel on the left, scalar 
channel on the right. 


CHARM 

Charmonium has been one of the most intensely studied probes in heavy-ion collisions since charmonium suppression 
was proposed as a signature of the QGP by Matsui and Satz [16]. There have been a number of lattice calculations of 
charmonium spectral functions, mostly in the quenched approximation [17, 18, 19, 20], but also recently using 2 h- 1 
light flavours [21]. We have previously studied charmonium on our 1st generation {Nf = 2) ensembles [22, 6]; here 
we will present results from our 2nd generation ensembles, as well as updated results for the momentum dependence 
of charmonium correlators from the 1st generation ensembles. 


Charmonium at zero and nonzero momentum 

Figure 3 shows charmonium spectral functions in the pseudoscalar (S-wave) and scalar (P-wave) channels from 
our second generation ensembles. Our results suggest that the S-wave ground state survives up to T « 1.67];, while 
the P-wave ground state melts close to Tc- This is consistent with our previous first generation results [22]. We have 
studied the dependence of these results on the default model used, and found that the results are stable for all but the 
highest temperature. We are in the process of analysing the data using the alternative Bayesian method of [2]. 

At nonzero momentum, the MEM reconstruction has much larger systematic uncertainties than at zero momentum 
[6, 23]. To circumvent this problem, we may instead compare the correlators at temperature T with the reconstructed 
correlators which result from integrating (1) using the kernel evaluated at T, but with a spectral function p{co,p-,Tr) 
determined at a reference temperature 7]. The reference temperature is chosen to be relatively low so that the MEM 
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FIGURE 4. Reconstructed charmonium correlators in the longitudinal (top) and transverse (bottom) vector channel (1st genera¬ 
tion), at two different momenta. The reference temperature is 7). = 230MeV (A^j- = 32). 


determination of p {co,p, Tr) is under control. Figure 4 shows this comparison for the vector channel at two different 
momenta with T) = 230MeV. Note that the comparison is always made with the reconstructed spectral function for the 
same momentum. We see that the combined effect of the temperature and momentum is different in the longitudinal 
and transverse channel: while the thermal modihcations in the longitudinal channel become smaller with increasing 
momentum, in the transverse channel they increase with the momentum, and are also larger in magnitude. 


Charmonium potential 

Nonrelativistic potential models have been widely used to study bound states of charm and beauty quarks both 
at zero and nonzero temperature. While it has long been known how their use can be given a hrm foundation at 
zero temperature in effective held theory (pNRQCD), it is only recently that a similar understanding has begun to be 
developed in the high-temperature case [24, 25]. Most if not all high-temperature applications have so far been based 
on the potential between inhnitely heavy (static) quarks, and it is not yet known how a hnite quark mass will modify 
this. 

We will take an alternative approach, assuming that the charmonium system can be described using a Schrodinger 
equation with a real potential, and determining the potential from the charmonium correlators that we compute on the 
lattice. Specihcally, we compute point-split correlators with an operator F = { 7 ^, 75 } 

Cr{r, t) = T)t/ (x,x + r)r\i/{x + r, t)v7(0, 0)r V(0, 0)). 
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FIGURE 5. Charmonium potential. Left: central (spin-independent) potential; right: spin-dependent potential. The points are 
horizontally shifted for clarity. The two sets of error bars (to the left and right of the symbols) denote statistical and systematic 
(from variations in the time range used) uncertainties, respectively. 




The potential yr(^) can then be determined using [26] 


dCr{r,z) 

dr 


\mc dP- 



Cr(r,T), 
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which is valid for T <C 1 jTT. The vector and pseudoscalar potentials can be combined to produce the central (spin- 
independent) and spin-dependent potentials Vc and Vs, 

Vcir) = (r) + (5), Vs{r) = Vy^ (r) - yy 3 (r). (6) 

The results are shown in hgure 5. At low temperatures the central potential is well represented by the conhning 
Cornell (linear + Coulomb) potential, while above Tc it becomes screened. The spin-dependent potential is attractive 
at intermediate distances with indications of a repulsive core. There appears to be a strong temperature dependence, 
also seen in our hrst generation results [27], which at present we have no clear understanding of. 


BEAUTY 

In recent years there has been an increased interest in b physics at high temperature and in heavy-ion collisions. In part 
this is because the energies reached at the LHC are such that bh pairs are created abundantly, whereas at RHIC and 
SPS they were relatively rare. Indeed, one of the early headline results from the heavy-ion collisions at CERN was the 
observation of sequential T suppression by the CMS collaboration [28]. 

Theoretically, the beauty system provides a cleaner probe of the QGP than the charm system, since effects such 
as regeneration and cold nuclear matter effects which may obscure the interpretation of the yields are much less 
significant. In addition, while the nonrelativistic approximation is marginal for charm quarks, it is clearly valid for 
bound states of beauty quarks. This allows us to use non-relativistic QCD (NRQCD), which is an effective theory 
obtained by integrating out the largest scale in the system, the heavy quark mass. In this case the kernel K in (1) 
simplihes to 

Kico,T) = e-‘^\ (7) 

Note that the heavy quark is explicitly not in thermal equilibrium here, and hence thermal (periodic) boundary 
conditions are not imposed. This has the additional advantage of doubling the number of independent points in the 
temporal direction. 

We have previously applied the NRQCD formalism to our hrst generation ensembles in a series of papers [29, 30, 
31, 32]. Here we will show results from our second generation ensembles [33]. The details of the simulation including 
the NRQCD action and the MEM reconstruction are presented in [33]. 

Our results [33] using the standard Bryan’s implementation of MEM are shown in hgure 6. We see that the S-wave 
(vector) ground state Y survives up to the highest temperatures studied, while there is no evidence for any surviving P- 
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FIGURE 6. Beautonium spectral functions in the vector channel (left) and axial-vector channel (right), using the maximum 
entropy method [33]. 




FIGURE 7. Beautonium spectral functions in the vector channel (left) and axial-vector channel (right), using the novel Bayesian 
method of [2]. 


wave (axial-vector) ground state for temperatures above Tc. This is consistent with our results from the first generation 
ensembles [30, 32]. 

Figure 7 shows our preliminary results using the novel Bayesian method of [2]. For the vector (S-wave) channel the 
results are in qualitative agreement with those from the maximum entropy method, although there are quantitative 
differences in the width of the ground state peak. For the axial-vector (P-wave) channel, on the other hand, the 
results in figure 7 suggest that the ground state survives until well into the QGP phase. This is in accordance with 
the results obtained in an analogous study using HotQCD ensembles [34]. We are currently investigating the source 
and signihcance of these differences, and this is discussed further in [35]. 


SUMMARY AND OUTLOOK 

We have presented results from the FASTSUM collaboration’s studies of high-temperature QCD using anisotropic 
lattices. Results have been obtained in the light quark, charm and beauty sectors, including the hrst results for the 
electrical conductivity and charge diffusion as a function of temperature below and above the phase transition, as well 
as charmonium and beautonium S and P wave states and the potential between charm quarks. 

In the beauty sector, we are currently working on a detailed comparison of the two different Bayesian methods 
used in the spectral reconstructions in order to fully understand the differences between them. We are also looking 
at alternative methods which can complement this understanding. We hope that this will lead to quantitative results 
for the temperature-dependent mass shift and width of T states, which will assist in the interpretation of results from 



































heavy-ion collisions. In the charm sector, we have obtained preliminary results for D mesons at high temperature, and 
we are also planning to study charm diffusion using the same methods as for the conductivity and charge diffusion 
(see [23] for preliminary results). We are also investigating baryons at high temperature, as well as the real-time static 
quark potential using the methods of [2, 36]. 

Finally, we are in the process of generating new ensembles with the same quark masses and spatial lattice spacing, 
but with a smaller temporal lattice spacing. The finer temporal resolution will lead to a more reliable spectral 
reconstruction and bring some of the main systematic uncertainties associated with this under control. It will also 
allow us to reach higher temperatures, and will be the first step towards the aim of providing quantitative predictions 
for spectral and transport properties of the QGP, that is results in the continuum limit with physical or near-physical 
quark masses. 
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